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0> ' Abstract 
<N . 

We present numerical evidence that a domain wall in a background with 
<^ , varying vacuum energy density acquires velocity in the direction of decreasing 

00 ' Hubble parameter. This should lead to at least a partial relaxation of the 

. cosmological constant on the wall. 

(N 
l> 
O , 

^ . Brane-world scenarios, according to which we live on a domain wall in a higher- 

dimensional space [|l|, open new avenues for thinking about the cosmological constant. 
Here, we want to explore feasibility of one particular mechanism that may help the 
Q_i' "constant" to relax to a small value. 

^ ■ The effective cosmological constant seen by an observer on the wall is determined 

by the local value of the Hubble parameter, H. It has long been known p[ that a 
domain wall in a space with zero vacuum energy will inflate, i.e. H will be nonzero. To 
^ I achieve H = 0, one may introduce a negative vacuum energy in the bulk, with a value 

; tuned precisely to compensate for the effect of the wall tension |0. Alternatively, one 

can introduce a scalar field with a dilaton-like coupling to the brane, to obtain static 
solutions for more generic values of the tension |Q . Or, with more than two external 
dimensions, one may keep the bulk vacuum energy zero but assume that gravity in 
the bulk is so weak that the Hubble parameter induced by the wall tension is small 
enough to be compatible with observations 

Here, we want to consider a different idea. If the bulk vacuum energy varies along 
an extra dimension, the wall may be moving in such a way that the Hubble parameter 
on it starts at a large positive value but then approaches zero — or even crosses zero 
and becomes negative. In this way, the wall will move towards a region where two 
large contributions to H (one from the bulk energy, and the other from the wall 
tension) nearly cancel. In other words, a finely tuned solution similar to that of ref. 
H will be reached automatically in the course of the evolution. If the cancellation 
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persists for some time, it may be possible to identify that period of time with the 
standard cosmology.^ 

How closely H can get to zero in this scenario will depend on the complicated 
dynamics occurring when the above two contributions become of the same order. 
In this paper, we address an earlier, simpler stage of the evolution, when the effect 
of the wall tension is still much smaller than the effect of the bulk vacuum energy. 
We have tried to obtain numerically geometries in which H changes little with time 
but varies significantly in space, and see which way a test domain wall will move 
in such a background. Our main result is that, in the backgrounds we were able 
to obtain, the wall quite generically acquires velocity in the direction of decreasing 
Hubble parameter, i.e. towards regions that are "less inflationary". Although the 
actual displacements that we have been able to detect in the lifetime of our simulation 
are quite small, it is encouraging that the final velocity is a sizable fraction of the 
speed of light. 

A suitable smooth background can be formed, for example, by another domain 
wall, which is much thicker than ours. The difference in the thickness may be due to a 
difference in masses, or simply in initial profiles, of the scalar fields forming the walls. 
Towards the end of the paper, we will speculate that there is also a more "economical" 
possibility, wherein a configuration of one domain wall in the background of another 
is realized with a single field in the context of topological inflation 0. For most of 
the paper, though, we will prefer to think about these two walls as being due to two 
different fields. 

In the remainder of the paper, we (i) describe some a priori guesses about geome- 
tries induced by a space-dependent vacuum energy and (ii) present a sample of the 
numerical results. 

We prefer to use isotropic coordinates, in which the metric has the form 

ds^ = -N\z,t)df + ^\z,t)[dz^ + dx^] ; (1) 

is called lapse, and ip — conformal factor. This choice of coordinates was adopted 
from our numerical studies of post-inflationary universe and the numerical al- 
gorithm we use here follows closely the one used in that work. We allow for the 
total of D > 1 spatial coordinates, but the metric depends only on one of them — the 
fifth coordinate z. Thus, rfx^ is the fiat Euclidean metric in the remaining [D — 1) 
dimensional space. 

For a preliminary discussion, it will suffice to use the energy constraint 

D-1 D-1 2 ^ 



* Similarly, H can change sign if the wall itself is static but the bulk vacuum energy or the wall 
tension change in time. In this case, however, it may be more difficult to make the time when H is 
close to zero long enough to include the entire known cosmological history. 
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and the equation for the lapse 



Here E is the energy density (defined below), and 

H = ^ (4) 

Nip ^ ' 

is the local Hubble parameter. The sum on the right-hand side of is over all scalar 
fields (pn of the model; a prime denotes derivative with respect to z, and a dot — with 
respect to t; = SnG; G is bulk Newton's constant. 

We first look for purely vacuum solutions with a constant energy density. In this 
case, the right-hand side of is zero, while E in equals the cosmological constant: 
E = A. The equations have the following solution: 

^ z — at 

N = ij\ (6) 

where A and a are constants. The Hubble parameter (|^) is constant: H = a/A^, 
so whenever a > the solution is infiationary. Note that for a > 0, the solution is 
well-defined in the half-space 2; > for any t < 0. This is appropriate, since according 
to (^, t is the conformal time, which for infiationary solutions is limited to negative 
values. 

Eq. (|) is satisfied by identically, but the energy constraint (0) imposes a 

relation between the parameters: 

This is the same relation as obtained for a class of solutions in ref. P], so ours 
should be the same solutions except written in different coordinates. In particular, 
for A = 0, we have a = ±1 and any A. [A will be determined by the junction 
condition on the wall: e.g. for a static thin wall with the geometry symmetrical 
about it, A~'^ = K,'^a/2{D — 1), where a is the wall tension.) The a = 1 solution is 
the infiationary solution of ref. p| in different coordinates.^] To obtain a = (and 
consequently H = 0), we need to have A = —D{D — 1)/2k^A^. This is the anti-de 
Sitter solution of ref. 0. 

Now, when A is not exactly a constant but is slowly varying, we expect the ge- 
ometry to interpolate between faster expanding regions with larger A and slower 

^Thc explicit relation between our coordinates and those used in ref. Q and called here z and i 
is: z — exp(— t) sinhz, t = — exp(— t) coshz, all in units of the constant Hubble distance H^^. 
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expanding, or perhaps even contracting, regions with smaller A. Geometries interpo- 
lating between de Sitter and and anti-de Sitter spaces can presumably be obtained by 
analytical continuation of certain Euclidean instantons. Such instantons, containing a 
de Sitter region sandwiched between two anti-de Sitter regions, have been considered 
in ref. in connection with the problem of counting de Sitter microstates. Note 
that for every positive a (expanding) solution, eq. (^ has a negative a (contract- 
ing) solution. Such contacting regions are ubiquitous in our numerical simulations. 
We do not know, however, if they can obtained by analytical continuation from any 
instantons. 

To study real-time evolution that may possibly lead to such interpolating geome- 
tries, we will need the rest of Einstein's equation. As in ref. 0, we write these in the 
Hamiltonian form. The Hamiltonian pair to eq. (^) is 

iV-"2"^ +2p-2)^ + ^ + ^^(2^-^ + ^^i^) • (8) 

(Although we did not emphasize that above, for a constant E = V = A eq. (§) is 
also satisfied by the solution (|)-(|^).) Equations for the fields are 

= (9) 



(^^^-^iV0;)'-iV^^^|^, (10) 



TTr, 



The energy density is 



V is the potential of the fields. The energy constraint Ce can be added to the right- 
hand side of d^) with any coefficient w (however, for numerical studies some choices 
may be better than others 0). In what follows, we use w = 1. 

The Hamiltonian structure of the equations allows us to use a leap-frog algorithm, 
in which the coordinates 0„ and are defined at full time steps, while the momenta 
TTn and H are defined at half-steps. The lapse does not have a Hamiltonian pair; as 
in 1^, we update it at full steps using eq. (^) and extrapolate it to half-steps. The 
energy constraint (^ is not used for the evolution but is monitored to provide an 
idea of how accurate the results are. In addition, we have monitored the momentum 



constraint 



n'^^J'CM = E ^n<Pn + iD- l)H' = (12) 

r n 



and have found that for the time intervals shown in the plots both Ce and Cm stay 
close to zero to a satisfactory degree. 
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To start numerical evolution, we need to solve the constraints (0) and ([T^) at the 



initial time slice. The momentum constraint (12) is solved trivially with H = const 
and the initial choice of the fields such that for each n either 7r„ or 0^ vanish. There 
is some advantage in being able to choose a suitable initial profile of ip by hand. To 
make our choice compatible with the energy constraint (^, we introduce a massive 
"radiation" field which is initially zero and whose initial momentum is determined 
by solving eq. (0). (This is similar to the method used in ref. [TO] to study the 



onset of infiation in inhomogeneous spacetimes.) We call this field "radiation" only 
to distinguish it from the fields forming the domain walls; its mass is in fact the 
largest mass in the problem. The "radiation" energy redshifts away in the course of 
the evolution (of course, more so in regions that are infiationary). 

We thus have three fields in our simulation: 0o that forms the thick background 
wall, 01 that forms the movable thin wall, and 02, the massive "radiation" field. We 
keep the thick wall centered at ^ = by choosing the boundary condition 0o(O, t) = 0. 
For the potential, we use 



^ = E ( -lf^l<Pl + \K<Pt^) + ^ + ^M^^ + A , (13) 



n=0,l 



where A is a constant — it is the vacuum energy at 2; = 0. A nontrivial z-dependence 
of vacuum energy is due to variation of the field 0o. With potential (p!3D, the only 
interaction between different fields is gravitational. 

A word about the choice of parameters. We measure times and lengths in units 
of /Xq^ and energy densities in units of fi^/Xo, so that we can set /io = 1 and Aq = 1. 
The results below are for /ii = 2, Ai = 2000, M = 250, = 37.5 (this dimensionless 

is k^/Iq/Ao in the original units), A = 0.16, and = 4 (i.e. one extra dimension). 

We want the background to be sufficiently smooth, so that the interior of the 
thick wall — a region near z = — can inflate. This type of inflation is sometimes 
called topological inflation So, we choose the initial proflle of 0o in the form 

V -2 + 

with the spatial scale i much larger than the wall's "natural" thickness: ip'^i ^ /Iq ^. 
The plots below are for i = L = 50, where L is the total length of the integration 
region: z G [0, L], from a grid with the total of 2049 points. 

Another necessary condition for topological inflation is that near z = the Hubble 
parameter supported by the vacuum energy of 0o is large enough, in comparison 
with fiQ 0. With the above values of the parameters, this Hubble parameter is 
Ho = (/t2 A/6) ^ I There have been numerical determinations of the critical ratios 
Hq/hq for which the interiors of various types of topological defects will inflate in 



three spatial dimensions |11|. We have not attempted a similarly detailed study for 



D = A, but have observed large amounts of inflation for the above value of Hq/ fiQ = 1. 
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Figure 1: Initial profiles of the two domain walls and of the conformal factor ijj. 

In comparison, gravitational effects due to the wall formed by 0i are, by our 
choice of the parameters, much smaller. A convenient measure of these effects is the 
expansion rate the wall would have if it were inflating in empty space with zero bulk 
cosmological constant: Hi = K^ai/2{D — 1), where o"i ~ yUi/Ai is the wall tension. 
For the above values of the parameters, this gives ai ~ 4 x 10^^. This estimate is 
well born out numerically,^ so we use Hi = 0.03. Because Hi <^ Hq, the wall formed 
by 01 has only a minor effect on the geometry in the rapidly expanding region near 
z = 0. On the other hand, when the wall is in a region where the expansion rate 
is of order Hi, or smaller, we expect strong gravitational effects due to the wall's 
tension. As we have already mentioned, such effects are likely to determine if the 
present scenario leads to a complete solution to the cosmological constant problem, 
but they lie outside the scope of this initial study. 

In Fig. |l], we show the initial profiles of (po, 0i, and ip for a typical initial position 
of the second wall, zi = 20. 

Now, let us see if an interpolating geometry of interest to us is reached in the 
course of numerical evolution. Evidence for that is shown in Fig. ^, where we plot 
profiles of the Hubble parameter H at several moments of time. We see that, after a 
transient behavior, H acquires a nearly static profile, in which it has positive values 
near z = (where the vacuum energy is the largest) and negative values at large z. 

To estimate the amount of inflation we obtain near 2; = 0, we plot in Fig. |^ the 
value of the conformal factor ip at z = a.s a function of time. We see an exponential 

■^-Variation of the tension with time is relatively small, which implies that the wall formed by (f>i 
maintains more or less constant physical (as opposed to comoving) thickness. 
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Figure 2: Profiles of tlie local Hubble parameter at several moments of time. 

growth by a factor over 40 between t = 12.5 and t = 20. This corresponds to the size 
of the universe growing by a factor over 1600.| 

We define velocity of the domain wall formed by 0i as the ratio of the momentum 
of 01 to its energy, both taken per unit transverse area: 



(15) 



where 



2 



^i-ai.+^+7Ai(0?-^i . (16) 



2v^ + ^ + y"' ~ a;, 

(The denominator in (^) divided by the 7-factor of the wall is our definition of the 
wall tension cxi.) Note that if the profile of 0i with respect to the comoving coordinate 
z does not change, i.e. dt(f)i = 0, the definition ([ISD gives zero. In other words, v as 
defined by (|I5|) is insensitive to a motion that is due solely to the expansion of the 
universe: it measures velocity of the wall relative to the expanding background. 

In Fig. ^ we plot the wall's velocity as a function of time, for the same initial 
position of the wall zi = 20 that was used for the previous plots, as well as for two 

§To properly interpret the plot, one needs to know also the time-dependence of the lapse N at 
z — 0, since it is N{z, t)dt that determines the increment of cosmic time at location z. As it happens, 
after iV reaches the value of 1 at z = in our simulation, we keep it equal to 1 there at all subsequent 
times (which is always possible due to the remaining gauge freedom). For Fig. this resetting of 
the clock rate took place at i « 12.7; after that, time t coincides (up to a constant) with the cosmic 
time at z — Q. 
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Figure 3: Inflation at the origin: conformal factor at 2; = as a function of time. 

smaller values of Zi. We see that v invariably becomes positive at late times, when a 
geometry with an almost static profile of the Hubble parameter is reached. The final 
value of the velocity can be a sizable fraction of the speed of light and is larger for 
walls that are further away from the origin. 

It may appear from Fig. ^ that the velocity goes to a steady asymptotic value, 
but we should point out that the lapse N{z, t) at the locations of the wall drops 
precipitously during the late-time evolution. So, the entire stage when v is close to 
its final value corresponds only to a small interval of cosmic time on the wall. As 
a consequence, the corresponding displacements of the wall are tiny. On the other 
hand, one can argue that the velocity remains steady over large amounts of expansion 
in the inflating region, cf. Fig. ^ 

Throughout this paper, we have considered the case when the wall formed by 0i 
is much "weaker", in terms of its gravitational effects, than the wall formed by 0o- 
This does not always have to be so, and, as we have already mentioned, perhaps the 
most interesting case with regard to the cosmological constant problem is when the 
"strengths" of the two walls are comparable. In addition, walls of different thicknesses 
but of equal "strengths" naturally arise in the context of topological inflation p with 
an extra dimension. Let us see how. 

The wall formed by 0i in our simulations maintains more or less constant physical 
thickness, as opposed to typical walls of 0o considered in topological inflation. Those 
latter are produced by growth of small fluctuations in the middle of the inflating 
region, where 0o ~ 0, and the vacuum energy E'vac is much larger than the gradient 
energy -Egrad of a new wall (whose thickness at birth is of the order of the Hubble dis- 
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Figure 4: Velocity of the domain wall formed by 0i as a function of time, for three 
different values of the initial position. 

tance). So, these walls expand in all directions 0. Analogy with the case considered 
here is much closer for a wall produced at the fringe of the inflating region. If a wall 
finds itself in a region with E^i^c ~ -Egrad, its thickness will not grow, but its gravi- 
tational effects will be comparable to those of the background. This can be seen by 
comparing the previously introduced parameters Hq and Hi. Indeed, Hq ~ k^E^^jc^ 
while Hi ~ K?ai ~ n'^EgradHQ^ , so Hi ~ Hq. 

Production of new domains walls of 0o should be observable even in the model 
considered here, if we were able to follow the evolution for a long enough time. In 
fact, it may well be that the eventual breakdown of our simulations (at times larger 
than those shown in the plots) is due primarily to the onset of this new, more com- 
plex dynamics. Obtaining good quality data for such a regime will require further 
numerical work. We hope to return to this question in a future publication. 

In conclusion, we have found that numerically solving Einstein's equations on one- 
dimensional lattices allows one to obtain, and experiment with, geometries in which 
the vacuum energy density varies along an extra dimension, and the local value of the 
Hubble parameter is almost static but varies from point to point. We have also found 
that a test domain wall in such a geometry will move in the direction of decreasing 
Hubble parameter, suggesting that this mechanism can lead to at least a partial 
relaxation of the cosmological "constant" in a brane world. 

This work was supported in part by the U.S. Department of Energy through Grant 
DE-FG02-91ER40681 (Task B). 
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